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ABSTRACT 

The relative efficiencies of the chemical pathways that can lead to the destruction of 
LiH and LiH"*" molecules, conjectured to be present in the primordial gas and to control 



^ . molecular cooling processes in the gravitational collapse of the post-recombination era, 

are revisited by using accurate quantum calculations for the several reactions involved. 
The new rates are employed to survey the behavior of the relative abundance of these 
I molecules at redshifts of interest for early universe conditions. We find significant differ- 

ences with respect to previous calculations, the present ones yielding LIH abundances 
■ higher than LiH"*" at all redshifts. 

O 

1. Introduction 



The suggestion that chemical processes involving lithium could play a role in the evolution of 
the early universe in the post-recombination era has been put forward several years ago (Lepp & 
'<^: Shull 1983; Standi et al. 1996, hereafter SLD96; Gahi & Palla 1998; hereafter GP98), starting from 

the consideration that Li is produced a few minutes after the Big Bang and the formation of other 
light atoms like H, D and He (Wagoner et al. 1967, Peebles 1993). The fractional abundances of 
these elements are sensitive to the baryon density of the universe and set constraints on its actual 
value (Cyburt et al. 2008). As the universe expanded, its radiation temperature decreased and 
the atomic ions originating from the above elements gave rise to neutral atoms by recombination 
with electrons, thus initiating the formation of molecular species by radiative association: H2, HD 
and LiH. The latter molecule, because of its large dipole moment and low ionization potential, may 
induce spatial and/or spectral distortions in the cosmic background radiation (CBR), as originally 
suggested by Dubrovich (1993) and as experimentally surmised by the pionering work of Maoli et 
al. (1994). Eventually, during the gravitational collapse leading to the formation of the first stars. 
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the low excitation threshold and the efficient radiative decay along the rotovibrational manifold 
are additional properties that could favor the role of LiH and LiH"'" as molecular coolants of the 
primordial gas (Bougleux & Galli 1997, hereafter BG97). The nonequilibrium level population 
in the presence of a different gas and radiation temperature may have a possible signature in 
protogalactic clouds, imprinting small fluctuations in the CBR spectrum (Schleicher et al. 2008), 
as we shall further discuss in our conclusions. 

The chemistry of Li in the early universe has been discussed in the past, reaching contrasting 
conclusions (see e.g., SLD96; BG97; GP98; Vonlanthen et al. 2009). Of critical relevance, is the 
uncertainty in the knowledge of reliable reaction rates for the destruction of LiH and LiH+ molecules 
via strongly exothermic reactions without entrance barriers: 

LiH + H Li + H2, (1) 

LiH + H+ ^ Li + H+, (2) 

and 

LiH+ + H ^ Li+ + Ha. (3) 

Therefore, it is an accurate knowledge of the reaction rates for the above processes, at low redshift 
values, that can ultimately tell us what the cnd-rolc of the LiH/LiH+ systems could be as efficient 
coolants under early universe conditions. The task of the present work is to show that the reaction 
rates recently determined from fully ab-initio quantum methods (Bovino et al. 2009, 2010a, 2010b), 
which also employ accurate interaction forces between partners, have a significant impact on the 
evolution of LiH and LiH+ during the post-recombination era of the early universe. We shall further 
show that a more realistic computational description of the rates for a neutralization process 

LiH+ + e~ ^> Li + H (4) 

could substantially change the relative abundance of the ionic molecular species. 

This paper is organized as follows: the details of the quantum methods for the new chemical 
reaction rates are given in Section 2, while Section 3 reports the list of considered reactions and 
describes the equations controlling the evolution of the gas temperature. In Section 4 we present 
the fractional abundances and an analysis of their behavior, while Section 5 summarizes our present 
conclusions. 



2. The quantum reactive calculations 

The study of the chemical evolution of the main species of the primordial gas requires accurate 
rate coefficients for the most important destruction reactions. In previous works, the Langevin 
approximation for barricrlcss reactions has been the main method used to obtain the rate coefficients 
for those reactions. In the last decade, with the developing of more realistic methods for reaction 
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processes, a great deal of work has been done to study reactions (1), (2) and (3). For instance, 
time-dependent (Bulut et al. 2008, Defazio ct al. 2005) and quasi-classical trajectory (Pino ct 
al. 2008, Dunne et al. 2001) calculations have been reported, and some of the calculated rates 
have been included in the abundance calculations of Lepp et. al. (2002). However, most of these 
calculations are restricted to a narrow range of energies or are limited to the high energy regime. 
Therefore, the main improvement of the present calculations has been the use of new, very accurate, 
potential energy surfaces (Martinazzo ct al. 2003a, 2003b, Wernli et al. 2009) combined with 
accurate quantum reactive calculations. The quantum methods range from the fully ab-initio 
coupled channels (Bovino et al. 2009) to the use of a negative imaginary potential approach (Bovino 
et al. 2010a, 2010b) for the ionic systems. The computed reaction probabilities and cross sections 
have been reported in previous papers and we shall not repeat here the details of the calculations. 
The final state-to-all cross section obtained from the probability (Pa-^all) expressed as 

where E is the collisional energy, ja the molecular rotational angular momentum, J the total 
angular momentum, and is the wave vector defined as 

kl = ^{E-Sa), (6) 

and Ea is the a-channel energy. Rate coefficients are computed by averaging the appropriate reactive 
cross sections over a Boltzmann distribution of velocities for the incoming atom: 

HTs) = ^^^ff"{rf/^' <7a-^all{E) expi-E/kBT^)E dE (7) 

{I^BJ-g) Jo 

where Tg is the gas temperature, ks is the Boltzmann constant and n the reduced mass of the 
system. The newly computed rate coefficients for reactions (1), (2) and (3) are reported in Bovino 
et al. (2009, 2010a, 2010b) and shown in Figure 1. The corresponding temperature-dependent fits 
are given in Table 1, along with all the reactions rates for the Li species employed in the present 
work. 



3. The chemical network and the evolutionary modelling 

The evolution of the pregalactic gas is considered in the framework of a Friedmann cosmological 
model and the primordial abundances of the main gas components are taken from the standard 
big bang nucleosynthesis results (Smith et al. 1993). The numerical values of the cosmological 
parameters used in the calculation are obtained from WMAP5 data (Komatsu et al. 2009) and are 
listed in Table 2. 
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In order to calculate the abundances of LiH and LiH+, a set of differential coupled chemical 
rate equations of the form: 

= kiormrijUk - kdestrii + ■■■ (8) 

has been solved. In eq. (8), /cform and /cdest are the formation and destruction reaction rates as 
listed in Table 1, and rii is the number density of the reactant species i. The evolution of the gas 
temperature Tg is governed by the equation (see e.g. GP98): 

^ = +3k^^^^- ^)Compton + (r - AUl], (9) 

where the first term represents the adiabatic cooling associated with the expansion of the Universe, 
R being the scale factor. The other two terms represent, respectively, the net transfer of energy 
from the CBR to the gas (per unit time and unit volume) via Compton scattering of CBR photons 
and electrons, 

,r-A)c,.p.,.=n.^^^I2^M^, (10) 

irieC 

and via excitation and de-excitation of molecular transition 

(r - A)moi = '^nk ^(xjCij - XjCji)hvij, (11) 

k i>j 

where Cij and Cji are the collisional excitation and de-excitation coefficients and Xi are the frac- 
tional level populations: for more details of this model see BG97 and GP98. The radiation temper- 
ature is Tr{z) = To{l + z), where Tq is the present-day CBR temperature. The chemical/thermal 
network is then completed by the rate equation for the redshift 

^ = I (12) 

where 

H{z) = HoV^ril + z)^ + ^m{l + zf + J^k(1 + + ^A. (13) 

Here, Hq is the Hubble constant and 0^, J^m, ^^K, ^^A are density parameters. The density n{z) of 
baryons at redshift z is 

n{z) = J^b"cr(l + zf, (14) 

with ricr = 3Hq / {STrGniB.) being the critical density. The calculations have been carried out from 
z = 10^ down to z = 1 and are discussed in the next Section. 



4. Results and discussion 



As discussed in the Introduction, one of the novelties of the present work is the use, for some 
of the most important chemical processes, of rates obtained from accurate quantum calculations 
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and not from qualitative estimates. It is therefore important for the reactions involving LiH/LiH+ 
to recall the behavior of the ab-initio reactive rates (in the temperature range of interest) that 
have been obtained via the quantum methods outlined in Section 2. The data reported in Figure 1 
indicate the depletion rates for the disappearance of LiH molecules by exothermically reacting 
with the surrounding atomic hydrogen (top panel). The temperature behavior of the ab-initio 
data is clearly very different from that assumed by SLD96. The new calculations of Bovino et 
al. (2009) indicate a variation over the same range of nearly three orders of magnitude. The same 
occurs for the destruction reaction of LiH by interaction with protons (middle panel): the ab-initio 
calculations (Bovino et al. 2010b) show that the rate varies by about one order of magnitude over 
the same range of temperatures. Finally, the quantum rates associated with the destruction of 
LiH"*" by the surrounding hydrogen gas (bottom panel) confirm the weak temperature dependence 
reported earlier (Lepp et al. 2002), although they turn out to be about a factor of 4 larger than the 
estimate of SLD96. Once the new quantities are employed within the evolutionary scheme outlined 
in Section 3, they yield a new set of production/destruction rates of LiH and LiH"'" as a function 
of redshift, as shown in Figure 2 and 3. 



4.1. LiH chemistry 

The varioTis curves shown in Figure 2 follow the notation given in Table 1 to the most relevant 
processes: those labelled as pi, . . . ,p4 (solid lines) refer to the different reactions leading to LiH 
production, while those labelled di,. . . ,d4 (dashed lines) refer to the destruction processes. A look 
at the figure allows one to make the following comments: 

(i) The production by radiative recombination from the Li(^P) state (pi) clearly dominates 
at high redshifts, while the production by recombination from the Li(^S') state (p2) becomes the 
dominant route from z ~ 800 down to z ~ 80. Then, the channel of recombination driven by 
electron attachment (^4) takes over. However, we must note that the latter rate is not known from 
quantum calculations, but only estimated by SLD96, together with the estimates of LiH formation 
from H- (SLD96). 

(a) At the highest redshifts, the destruction channel is dominated by photodissociation (di), 
a rate estimated from the detailed balance of process (^2), while at lower z values {z < 300) the 
collisional process of reactive destruction discussed earlier (^2) takes over and dominates the whole 
destruction channel down to the lowest redshifts. The same reactive process with protons (1^3) is 
seen to be less important, although larger than previous estimates (^4). 

In conclusion, the above data indicate that the radiative photoassociation from the n = 25 level 
of Li is the most significant formation path, while the collisional destruction by H (^2) overcomes 
production over the whole range of redshifts below 300. 
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4.2. LiH+ chemistry 

Figure 3 shows the processes leading to the formation and the destruction of LiH+, labelled as 
p'l,. . . ,P2 and d^, . . . , ^4" as in Table 1. As clearly shown by the figure, the radiative association 
processes (p^ and P2) are seen to dominate the formation of LiH"*" over a very broad range of 
redshifts. The association involving the Li ion (p^) is the major pathway down to about z c:^ 700, 
when the reaction of with neutral Li (p^) takes over and becomes dominant by several orders 
of magnitude over the rates of the charge-exchange reaction (p^) estimated by SLD96. 

As for the destruction processes for LiH"^, one sees that photodissociation (df) is the dominant 
destruction path of the ionic molecule, due to its small dissociation energy value. Furthermore, the 
dominant dissociation channel comes from Li"*" production (d^), while the one yielding H"*" (^2^) is 
relevant only at the highest redshifts. The electron-induced dissociation (^4") was earlier based on a 
simple estimates from SLD96 while now we have included the more accurate quantum calculations 
from Curik Sz Greene (2007,2008) which yield much larger rates. It turns out to be important only 
at low redshifts. On the other hand, the chemical paths which involve the reactions that destroy 
LiH+ (d^) are seen to be indeed the most important processes around 2; ~ 40 while becoming 
less efficient than the d'^ process below z ~ 10. The electron-assisted dissociation of LiH+ {d^) is 
therefore important only at the lowest redshifts. 

In summary, the ionic formation is dominated by the radiative association of Li with H"*" over a 
very broad range of redshifts, while destruction by photodissociation prevails at larger redshift and 
chemical destruction (dg") dominates at lower redshifts. However, down to 2; ~ 1 the dissociative 
neutralization (d^) is markedly dominant. 

4.3. Chemical evolution 

The new evolution of relative abundances of atomic and molecular species involving Li is 
displayed in Figure 4, comparing the earlier results obtained with the chemical network of GP98 
(dashed lines) and the present results (solid lines) which employ the new quantum rates. The main 
results are the following: 

(i) The temperature dependence of the quantum rates, shown by the data of Figure 1, is clearly 
reflected in the molecular abundances of LiH at redshifts below z ~ 300. The rapid rise of LiH at 
z ~ 300 found by GP98 is no longer seen owing to the increased destruction rate of channel (1). 
The rise at 2; ~ 100 is due to the importance of LiH formation by Li~ that persists at all redshifts 
together with the radiative association of Li with H. This is drastically different from the behavior 
found by GP98 where the constant value of the rate for reaction (1) caused the steady drop of LiH. 
The final {z = 0) value of LiH is ~ 7 x 10"^*^ with an increase of a factor of ~ 70 relative to GP98. 

(ii) Due to the differences in behavior shown by the cross sections in Figure 1, the ionic partner, 
LiH+, behaves differently at 2 < 30. First, the increased efficiency of LiH+ destruction (reaction 3; 
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rate d^) limits the sharp rise in abundance at z ~ 20. Then, the gentle decline below z ~ 5 is due 
to the efficiency of electronic recombination (rate df). The new final abundances of LiH"'" (solid 
line) is now smaller than earlier estimates of GP98 by a factor of ~ 20. 

The scenario emerging from the above calculations therefore indicates that in the regions of 
redshift below z ~ 30, LiH remains the more abundant species compared to LiH+, but only by a 
factor of ~ 2-7. The two molecular abundances, on the other hand, reach the largest values for 
z < 10, remaining both fairly small (LiH~10~^^, LiH+ ~10~^^) and hard to detect. 

5. Conclusions 

Building upon recent quantum reactive calculations (Bovino et al. 2009, 2010a, 2010b) in- 
volving the chemical evolution of LiH/LiH+ molecules within the expected conditions in the early 
universe, the present work has revisited the analysis of all the dynamical processes that are known 
to significantly contribute to the production/destruction of the lithium-containing molecules. We 

have employed as many results as possible from calculations based on ab-initio methods, both 
for the interaction forces and the quantum dynamics, resorting to estimates only for a few of the 
considered processes. 

One of the main results from the quantum reactive calculations is that two of the important 
rates exhibit a temperature dependence that was not present in the earlier estimates. 

The results for the abundance of LiH indicate that this molecule is now much more likely to 
have survived at low redshift, since its fractional abundance, albeit still fairly small, goes up by a 
factor of ^ 70 compared to previous estimates. However, it becomes smaller by nearly the same 
amount for z > 300. 

The fractional abundance of LiH+ only becomes significant in the low-redshift region of z < 30 
and shows a reduced value of the relative abundance by about one order of magnitude with respect 
to the earlier estimates. The comparison between the specific fractional abundances of the two 
species now indicates that the neutral molecule is likely to be more abundant than the ionic species 
and that their ratio in the region of small redshift increases up to a factor of 7 down to z ~ 1. Unlike 
previous estimates that predicted a difference at low-z of about two orders of magnitude in favor 
of LiH+, the new calculations find more comparable abundances of LiH and LiH+. Furthermore, 
we find that the neutral molecule, in spite of a large dilution at the low redshifts, should be more 
amenable to experimental observation (e.g., Persson et al. 2010) than its ionic counterpart. 

Finally, given the recent studies (e.g. see Zaldagarriaga and Loeb, 2002) on the possible detec- 
tion, wihin the microwave background anisotropics, of the earlier imprint from the recombination 
history of lithium, the present results strongly suggest that such anisotropics could be amenable to 
observation due to the changed optical depth induced by the changes on lithium abundances 
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Table 1: List of considered reactions. Labels in the second column refer to Figures 2 and 3. 

reaction rate (cm"^ or s""*^) notes reference 



1) 

2) 
3) 

4) 
S) 
6) 
7) 
8) 
9) 



10) 
11) 



12) 

13) 

14) 



15) 
16) 

17) 
18) 



PI 
PI 
P2 



d2 



d4 



19) d. 

20) d 

21) d. 

22) d 

23) d 

24) d 
26) d. 
26) 

27) 
28) 



l,irp) + H -> LiH + hv 
Li(2p) + H -> LiH + hf 
Li + H LiH + hv 



H 



+ 



LiH + o 
LiH+ + h 
+ H+ -> Li+ + H 
+ H+ -> Li+ + H ■ 
+ c ~> Li^ + hiy 
+ H± LiH + H+ 



Li+ 



» LiH+ + hv 
Li + hv 



H -^ LiH 
H+ Li 



LiH + H -) 
LiH + H+ 



Li + Ha 



LiH 
LiH 



H+ LiH+ + H 
H+ -> Li+ + H2 



LiH+ ■ 
LiH+ ■ 
LiH+ ■ 
LiH+ ■ 



■+ Li+ + H2 
■+ Li + H+ 
■* LiH + H+ 
. Li + H 



LiH + hi/ ^ Li + H 
LiH+ + ^ Li+ + H 
LiH+ + hi/ ^- Li + H+ 
Li + hi/ ^ Li+ + c 
Li~ + hi/ — f Li + e 
LiH + hi/ ->. Li(2p) + H 



Li reactions 
2.0 X 10~'^'^T°'l* exp(-Tg/B100) 
1.9 X io-llT-0.34 

4.0 X 10-2" ^ 

I.J -p J-JAoJ^"^^ 

fV-io 

-0.49 



exp[-Tg/40(i,T + (Tg/13193)^] 
4.0 X i(r "■ 
5.3 X lO^^TJ 



2.5 X 10 



'■9 exp(-Tg/121()) 



1.7 X io-13t~'' ''51 cxp{-Tg/282000) 

6.1 X IO-i^tI-^** oxp(-Tg/17200) 

6.3 X 10~^° exp(-2553/Tg) 

7.2 X 10~'^*T^-1* exp(-1470/Tg) 

Li^" reactions 

1.4 X 10-20^-0.9 exp{-Tg/7000) 

1.036 X 10-^*[-\/Tg/107.7x 
(1 + yTg/107.7)° '"'12 X 

(1 + yjr^/i.nr x io7)i 3388]-i 

6.3 X 10-''Tg"^'^2(j^ _,_ Tg/14000) 

Li" reactions 

4.0 X 10-^° 

2.3 X 10-''Tg~^/2 

LiH reactions 
2.0 X lQ^''-^Tg cicp(-Ti^/12t)0) 
2.9 X i()-lf)Tf)-"'0 
-2.6 X lO-lOTgO-SOexp 
1.0 X 10-^ 
1.0 X 10-^ 



i(-400/Tg) 



8. 


,7 


X 


10 


9. 


,0 


X 


10 


1. 


,0 


X 


10 


3. 


,9 


X 


10 



LiH+ reactions 
0^0.040 exp(Tg/5.92 X 10^) 
exp(-66400/Tg) 
^ exp(-67900/Tg) 
T-^ '^^ exp(-Tg/1200) 

photoreactions 



A^S+ -J- quantal calc. (a) 

B^n Xis+, quantal calc. (a) 

quantal calc. {b),(c),(d) 

estimate (e) 

quantal calc. (c),(d) 

quantal calc. (f) 

quantal calc. (g) 
det. bal. applied to (27) 

Tg < 500, quantal calc. (i) 

Tg > 500, quantal calc. (j) 



quantal calc. (c),(d) 

quantal calc. (k) 

Landau-Zener approx. (1) 

estimate (e) 

Landau-Zener approx. (1) 

quantal calc. (m) 

quantal calc. (n) 

quantal calc. (o) 

estimate (e) 



quantal calc. (p) 

estimate (e) 

estimate (e) 

quantal calc. (q) 



det. bal. applied to (3) 

det. bal. applied to (10) 

det. bal. applied to (5) 

det. bal. applied to (11) 

quantal calc. (h) 

det. bal. applied to (1), (2) 



(a) Gianturco Sz Gori Giorgi (1996); (b) Bennett et al. (2003, 2008); (c) Dalgarno, Kirby &z Stancil (1996); (d) Gianturco &z Gori Giorgi (1997); 
(e) Stancil, Lcpp &z Dalgarno (1996); (f) Kimura, Dutta ^ Shimakura (1994); (g) Stancil & Zygclman (1996); (h) Ramsbottom, Bell & 
Berrington (1994); (i) Bulut ct al. (2009); (j) Pino ct al. (2008); (k) Vcrncr &z Fcrland (1996); (1) Croft, Dickinson & Gadca (1999): (m) Bovino, 
Wernli & Gianturco (2009); (n) Bovino et al. (2010); (o) Bulut et al. (2008); (p) Bovino, Stoecklin & Gianturco (2010); (q) Curik & 
Greene (2007, 2008). 
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Fig. 1. — A comparison of reaction rates for: destruction of LiH with H reaction (top panel) and H"*" 
(middle panel) ; destruction of LiH"*" with H (lower panel) . The solid lines are fits of the computed 
data (open triangles) which used quantum reactive scattering methods (see Section 2). Dashed 
curved represent estimates by SLD96. 
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Fig. 2. — Computed production/destruction rates of LiH as function of redshift. Solid and dashed 
curves represent production and destruction processes, respectively. See text for details. 




Fig. 3. — Same as Figure 2 for LiH+. See text for details. 
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Fig. 4. — Relative abundances of Li, Li~, LiH, and LiH+ in the post-recombination era as function 

of redshift z (lower scale) and gast temperature Tg (upper scale): present results (solid curves), 
results of GP98 (dashed curves). The dotted lines at z=ll and z = 7 indicate the approximate 
redshifts at which reionization of the primordial gas starts {z ~ 11) and is completed {z ~ 7), 
according to Spergel et al. (2007). 



Table 2: Cosmological model 
Parameter Numerical value 



Ho 


70.5 km s"^ Mpc"^ 


To 


2.725 K 




3141 




0.228 




0.0456 




^dm ~l~ 




^m/ (1 "I" •^eq) 




0.726 


Qk 




/h 


0.924 


/He 


0.076 


/d 


2.38 X 10"^ 


/Li 


4.04 X 10"^° 



